Population genomic analyses reveal that salinity and geographic isolation drive diversification in a free-living protist.

Protists make up the vast diversity of eukaryotic life and play a critical role in biogeochemical cycling and in food webs. Because of their small size, cryptic life cycles, and large population sizes, our understanding of speciation in these organisms is very limited. We performed population genomic analyses on 153 strains isolated from eight populations of the recently radiated dinoflagellate genus Apocalathium, to explore the drivers and mechanisms of speciation processes. Species of this genus inhabit both freshwater and saline habitats, lakes and seas, and are found in cold temperate environments across the world. RAD sequencing analyses revealed that the populations were overall highly differentiated, but morphological similarity was not congruent with genetic similarity. While geographic isolation was to some extent coupled to genetic distance, this pattern was not consistent. Instead, we found evidence that the environment, specifically salinity, is a major factor in driving ecological speciation in Apocalathium. While saline populations were unique in loci coupled to genes involved in osmoregulation, freshwater populations appear to lack these. Our study highlights that adaptation to freshwater through loss of osmoregulatory genes may be an important speciation mechanism in free-living aquatic protists.

Local adaptation through countergradient selection in northern populations of Skeletonema marinoi.

Marine microorganisms have the potential to disperse widely with few obvious barriers to gene flow. However, among microalgae, several studies have demonstrated that species can be highly genetically structured with limited gene flow among populations, despite hydrographic connectivity. Ecological differentiation and local adaptation have been suggested as drivers of such population structure. Here we tested whether multiple strains from two genetically distinct Baltic Sea populations of the diatom Skeletonema marinoi showed evidence of local adaptation to their local environments: the estuarine Bothnian Sea and the marine Kattegat Sea. We performed reciprocal transplants of multiple strains between culture media based on water from the respective environments, and we also allowed competition between strains of estuarine and marine origin in both salinities. When grown alone, both marine and estuarine strains performed best in the high-salinity environment, and estuarine strains always grew faster than marine strains. This result suggests local adaptation through countergradient selection, that is, genetic effects counteract environmental effects. However, the higher growth rate of the estuarine strains appears to have a cost in the marine environment and when strains were allowed to compete, marine strains performed better than estuarine strains in the marine environment. Thus, other traits are likely to also affect fitness. We provide evidence that tolerance to pH could be involved and that estuarine strains that are adapted to a more fluctuating pH continue growing at higher pH than marine strains.

The new chytridiomycete Paradinomyces triforaminorum gen. et sp. nov. co-occurs with other parasitoids during a Kryptoperidinium foliaceum (Dinophyceae) bloom in the Baltic Sea.

A new chytrid genus and species was isolated and cultured from samples obtained in the Baltic Sea during a dinoflagellate bloom event. This species is characterized by having a spherical sporangium without papillae and zoospores of 2-3 µm in diameter that are released through 3 discharge pores. Molecular phylogeny based on ribosomal operon showed its sister position to the Dinomyces cluster in Rhizophydiales. Zoospores lack fenestrated cisternae but contain a paracrystalline inclusion, found in a Rhizophydiales representative for the first time. Additionally, the kinetid features are uncommon for Rhizophydiales and only observed in Dinomyces representatives so far. These morphological features and its phylogenetic relationships justify the description of the new genus and speciesParadinomyces triforaminorum gen. nov. sp. nov. belonging to the family Dinomycetaceae. The chytrid was detected during a high-biomass bloom of the dinoflagellate Kryptoperidinium foliaceum. Laboratory experiments suggest this species is highly specific and demonstrate the impact it can have on HAB development. The chytrid co-occurred with three other parasites belonging to Chytridiomycota (Fungi) and Perkinsea (Alveolata), highlighting that parasitic interactions are common during HABs in brackish and marine systems, and these multiple parasites compete for similar hosts.

Strain-specific transcriptional responses overshadow salinity effects in a marine diatom sampled along the Baltic Sea salinity cline.

The salinity gradient separating marine and freshwater environments represents a major ecological divide for microbiota, yet the mechanisms by which marine microbes have adapted to and ultimately diversified in freshwater environments are poorly understood. Here, we take advantage of a natural evolutionary experiment: the colonization of the brackish Baltic Sea by the ancestrally marine diatom Skeletonema marinoi. To understand how diatoms respond to low salinity, we characterized transcriptomic responses of acclimated S. marinoi grown in a common garden. Our experiment included eight strains from source populations spanning the Baltic Sea salinity cline. Gene expression analysis revealed that low salinities induced changes in the cellular metabolism of S. marinoi, including upregulation of photosynthesis and storage compound biosynthesis, increased nutrient demand, and a complex response to oxidative stress. However, the strain effect overshadowed the salinity effect, as strains differed significantly in their response, both regarding the strength and the strategy (direction of gene expression) of their response. The high degree of intraspecific variation in gene expression observed here highlights an important but often overlooked source of biological variation associated with how diatoms respond to environmental change.

Seasonal genotype dynamics of a marine dinoflagellate: Pelagic populations are homogeneous and as diverse as benthic seed banks.

Genetic diversity is the basis for evolutionary adaptation and selection under changing environmental conditions. Phytoplankton populations are genotypically diverse, can become genetically differentiated within small spatiotemporal scales and many species form resting stages. Resting stage accumulations in sediments (seed banks) are expected to serve as reservoirs for genetic information, but so far their role in maintaining phytoplankton diversity and in evolution has remained unclear. In this study we used the toxic dinoflagellate Alexandrium ostenfeldii (Dinophyceae) as a model organism to investigate if (i) the benthic seed bank is more diverse than the pelagic population and (ii) the pelagic population is seasonally differentiated. Resting stages (benthic) and plankton (pelagic) samples were collected at a coastal bloom site in the Baltic Sea, followed by cell isolation and genotyping using microsatellite markers (MS) and restriction site associated DNA sequencing (RAD). High clonal diversity (98%-100%) combined with intermediate to low gene diversity (0.58-0.03, depending on the marker) was found. Surprisingly, the benthic and pelagic fractions of the population were equally diverse, and the pelagic fraction was temporally homogeneous, despite seasonal fluctuation of environmental selection pressures. The results of this study suggest that continuous benthic-pelagic coupling, combined with frequent sexual reproduction, as indicated by persistent linkage equilibrium, prevent the dominance of single clonal lineages in a dynamic environment. Both processes harmonize the pelagic with the benthic population and thus prevent seasonal population differentiation. At the same time, frequent sexual reproduction and benthic-pelagic coupling maintain high clonal diversity in both habitats.

Cyanobacterial Akinete Distribution, Viability, and Cyanotoxin Records in Sediment Archives From the Northern Baltic Sea.

Cyanobacteria of the order Nostocales, including Baltic Sea bloom-forming taxa Nodularia spumigena, Aphanizomenon flosaquae, and Dolichospermum spp., produce resting stages, known as akinetes, under unfavorable conditions. These akinetes can persist in the sediment and germinate if favorable conditions return, simultaneously representing past blooms and possibly contributing to future bloom formation. The present study characterized cyanobacterial akinete survival, germination, and potential cyanotoxin production in brackish water sediment archives from coastal and open Gulf of Finland in order to understand recent bloom expansion, akinete persistence, and cyanobacteria life cycles in the northern Baltic Sea. Results showed that cyanobacterial akinetes can persist in and germinate from Northern Baltic Sea sediment up to >40 and >400 years old, at coastal and open-sea locations, respectively. Akinete abundance and viability decreased with age and depth of vertical sediment layers. The detection of potential microcystin and nodularin production from akinetes was minimal and restricted to the surface sediment layers. Phylogenetic analysis of culturable cyanobacteria from the coastal sediment core indicated that most strains likely belonged to the benthic genus Anabaena. Potentially planktonic species of Dolichospermum could only be revived from the near-surface layers of the sediment, corresponding to an estimated age of 1-3 years. Results of germination experiments supported the notion that akinetes do not play an equally significant role in the life cycles of all bloom-forming cyanobacteria in the Baltic Sea. Overall, there was minimal congruence between akinete abundance, cyanotoxin concentration, and the presence of cyanotoxin biosynthetic genes in either sediment core. Further research is recommended to accurately detect and quantify akinetes and cyanotoxin genes from brackish water sediment samples in order to further describe species-specific benthic archives of cyanobacteria.

Description of two new coexisting parasitoids of blooming dinoflagellates in the Baltic sea: Parvilucifera catillosa sp. nov. and Parvilucifera sp. (Perkinsea, Alveolata).

Perkinsea are a group of intracellular protist parasites that inhabit all types of aquatic environments and cause significant population declines of a wide variety of hosts. However, the diversity of this lineage is mostly represented by environmental rDNA sequences. Complete descriptions of Perkinsea that infect marine dinoflagellates have increased in recent literature due to the identification, isolation and culturing of representatives during bloom events, contributing to expand the knowledge on the diversity and ecology of the group. Shallow coastal areas in the Baltic Sea suffer seasonal dinoflagellate blooms. In summer 2016, two parasitoids were isolated during a Kryptoperidinium foliaceum bloom in the Baltic Sea. Morphological features and sequences of the small and large subunit of the ribosomal DNA gene revealed these two parasitoids were new species that belong to the genus Parvilucifera. This is the first time that Parvilucifera infections are reported in the Inner Baltic Sea. The first species, Parvilucifera sp. has some morphological and phylogenetic features in common with P. sinerae and P. corolla, although its ultrastructure could not be studied and the formal description could not be done. The second new species, named Parvilucifera catillosa, has several distinct morphological features in its zoospores (e.g. the presence of a rostrum), and in the shape and size of the apertures in the sporangium stage, which are larger and more protuberant than in the other species of the genus. Infections observed in the field and cross-infection experiments determined that the host range of both Parvilucifera species was restricted to dinoflagellates, each one showing a different host preference. The coexistence in the same environment by the two closely related parasitoids with very similar life cycles suggests that their niche separation is the preferred host.

Morphological and phylogenetic data do not support the split of Alexandrium into four genera.

A recently published study analyzed the phylogenetic relationship between the genera Centrodinium and Alexandrium, confirming an earlier publication showing the genus Alexandrium as paraphyletic. This most recent manuscript retained the genus Alexandrium, introduced a new genus Episemicolon, resurrected two genera, Gessnerium and Protogonyaulax, and stated that: "The polyphyly [sic] of Alexandrium is solved with the split into four genera". However, these reintroduced taxa were not based on monophyletic groups. Therefore this work, if accepted, would result in replacing a single paraphyletic taxon with several non-monophyletic ones. The morphological data presented for genus characterization also do not convincingly support taxa delimitations. The combination of weak molecular phylogenetics and the lack of diagnostic traits (i.e., autapomorphies) render the applicability of the concept of limited use. The proposal to split the genus Alexandrium on the basis of our current knowledge is rejected herein. The aim here is not to present an alternative analysis and revision, but to maintain Alexandrium. A better constructed and more phylogenetically accurate revision can and should wait until more complete evidence becomes available and there is a strong reason to revise the genus Alexandrium. The reasons are explained in detail by a review of the available molecular and morphological data for species of the genera Alexandrium and Centrodinium. In addition, cyst morphology and chemotaxonomy are discussed, and the need for integrative taxonomy is highlighted.

Dead or alive: sediment DNA archives as tools for tracking aquatic evolution and adaptation.

DNA can be preserved in marine and freshwater sediments both in bulk sediment and in intact, viable resting stages. Here, we assess the potential for combined use of ancient, environmental, DNA and timeseries of resurrected long-term dormant organisms, to reconstruct trophic interactions and evolutionary adaptation to changing environments. These new methods, coupled with independent evidence of biotic and abiotic forcing factors, can provide a holistic view of past ecosystems beyond that offered by standard palaeoecology, help us assess implications of ecological and molecular change for contemporary ecosystem functioning and services, and improve our ability to predict adaptation to environmental stress.

Generalist Life Cycle Aids Persistence of Alexandrium ostenfeldii (Dinophyceae) in Seasonal Coastal Habitats of the Baltic Sea.

In seasonal environments, strong gradients of environmental parameters can shape life cycles of phytoplankton. Depending on the rate of environmental fluctuation, specialist or generalist strategies may be favored, potentially affecting life cycle transitions. The present study examined life cycle transitions of the toxin producing Baltic dinoflagellate Alexandrium ostenfeldii and their regulation by environmental factors (temperature and nutrients). This investigation aimed to determine whether genetic recombination of different strains is required for resting cyst formation and whether newly formed cysts are dormant. Field data (temperature and salinity) and sediment surface samples were collected from a site with recurrent blooms and germination and encystment experiments were conducted under controlled laboratory conditions. Results indicate a lack of seasonal germination pattern, set by an endogenous rhythm, as commonly found with other dinoflagellates from the Baltic Sea. Germination of quiescent cysts was triggered by temperatures exceeding 10°C and combined nutrient limitation of nitrogen and phosphorus or a drop in temperature from 16 to 10°C triggered encystment most efficiently. Genetic recombination was not mandatory for the formation of resting cysts, but supported higher numbers of resistant cysts and enhanced germination capacity after a resting period. Findings from this study confirm that A. ostenfeldii follows a generalist germination and cyst formation strategy, driven by strong seasonality, which may support its persistence and possibly expansion in marginal environments in the future, if higher temperatures facilitate a longer growth season.

Exploring Molecular Signs of Sex in the Marine Diatom Skeletonema marinoi.

Sexual reproduction plays a fundamental role in diatom life cycles. It contributes to increasing genetic diversity through meiotic recombination and also represents the phase where large-sized cells are produced to counteract the cell size reduction process that characterizes these microalgae. With the aim to identify genes linked to the sexual phase of the centric planktonic diatom Skeletonemamarinoi, we carried out an RNA-seq experiment comparing the expression level of transcripts in sexualized cells with that of large cells not competent for sex. A set of genes involved in meiosis were found upregulated. Despite the fact that flagellate gametes were observed in the sample, we did not detect the expression of genes involved in the synthesis of flagella that were upregulated during sexual reproduction in another centric diatom. A comparison with the set of genes changing during the first phases of sexual reproduction of the pennate diatom Pseudo-nitzschiamultistriata revealed the existence of commonalities, including the strong upregulation of genes with an unknown function that we named Sex Induced Genes (SIG). Our results further broadened the panel of genes that can be used as a marker for sexual reproduction of diatoms, crucial for the interpretation of metatranscriptomic datasets.

Distributions of three Alexandrium species and their toxins across a salinity gradient suggest an increasing impact of GDA producing A. pseudogonyaulax in shallow brackish waters of Northern Europe.

Blooms of Alexandrium spp. are a well-known phenomenon in Northern European waters. While A. tamarense/catenella, and A. pseudogonyaulax have been reported from marine waters, high densities of A. ostenfeldii are mainly observed at lower salinities in North Sea estuaries and the Baltic Sea, suggesting salinity as a driver of Alexandrium species composition and toxin distribution. To investigate this relationship, an oceanographic expedition through a natural salinity gradient was conducted in June 2016 along the coasts of Denmark. Besides hydrographic data, phytoplankton and sediment samples were collected for analyses of Alexandrium spp. cell and cyst abundances, for toxin measurement and cell isolation. Plankton data revealed the predominance of A. pseudogonyaulax at all transect stations while A. ostenfeldii and A. catenella generally contributed a minor fraction to the Alexandrium community. High abundances of A. pseudogonyaulax in the shallow enclosed Limfjord were accompanied by high amounts of goniodomin A (GDA). This toxin was also detected at low abundances along with A. pseudogonyaulax in the North Sea and the Kattegat. Genetic and morphological characterization of established strains showed high similarity of the Northern European population to distant geographic populations. Despite low cell abundances of A. ostenfeldii, different profiles of cycloimines were measured in the North Sea and in the Limfjord. This field survey revealed that salinity alone does not determine Alexandrium species and toxin distribution, but emphasizes the importance of habitat conditions such as proximity to seed banks, shelter, and high nutrient concentrations. The results show that A. pseudogonyaulax has become a prominent member of the Alexandrium spp. community over the past decade in the study area. Analyses of long term monitoring data from the Limfjord confirmed a recent shift to A. pseudogonyaulax dominance. Cyst and toxin records of the species in Kiel Bight suggest a spreading potential into the brackish Baltic Sea, which might lead to an expansion of blooms under future climate conditions.

Genetic relatedness of a new Japanese isolates of Alexandrium ostenfeldii bloom population with global isolates.

In recent years, blooms of toxic Alexandrium ostenfeldii strains have been reported from around the world. In 2013, the species formed a red tide in a shallow lagoon in western Japan, which was the first report of the species in the area. To investigate the genetic relatedness of Japanese A. ostenfeldii and global isolates, the full-length SSU, ITS and LSU sequences were determined, and phylogenetic analyses were conducted for isolates from western and northern Japan and from the Baltic Sea. Genotyping and microsatellite sequence comparison were performed to estimate the divergence and connectivity between the populations from western Japan and the Baltic Sea. In all phylogenetic analyses, the isolates from western Japan grouped together with global isolates from shallow and low saline areas, such as the Baltic Sea, estuaries on the east coast of U.S.A. and from the Bohai Sea, China. In contrast, the isolates from northern Japan formed a well-supported separate group in the ITS and LSU phylogenies, indicating differentiation between the Japanese populations. This was further supported by the notable differentiation between the sequences of western and northern Japanese isolates, whereas the lowest differentiation was found between the western Japanese and Chinese isolates. Microsatellite genotyping revealed low genetic diversity in the western Japanese population, possibly explained by a recent introduction to the lagoon from where it was detected. The red tide recorded in the shallow lagoon followed notable changes in the salinity of the waterbody and phytoplankton composition, potentially facilitating the bloom of A. ostenfeldii.

Future temperature and salinity do not exert selection pressure on cyst germination of a toxic phytoplankton species.

Environmental conditions regulate the germination of phytoplankton resting stages. While some factors lead to synchronous germination, others stimulate germination of only a small fraction of the resting stages. This suggests that habitat filters may act on the germination level and thus affect selection of blooming strains. Benthic "seed banks" of the toxic dinoflagellate Alexandrium ostenfeldii from the Baltic Sea are genetically and phenotypically diverse, indicating a high potential for adaptation by selection on standing genetic variation. Here, we experimentally tested the role of climate-related salinity and temperature as selection filters during germination and subsequent establishment of A. ostenfeldii strains. A representative resting cyst population was isolated from sediment samples, and germination and reciprocal transplantation experiments were carried out, including four treatments: Average present day germination conditions and three potential future conditions: high temperature, low salinity, and high temperature in combination with low salinity. We found that the final germination success of A. ostenfeldii resting cysts was unaffected by temperature and salinity in the range tested. A high germination success of more than 80% in all treatments indicates that strains are not selected by temperature and salinity during germination, but selection becomes more important shortly after germination, in the vegetative stage of the life cycle. Moreover, strains were not adapted to germination conditions. Instead, highly plastic responses occurred after transplantation and significantly higher growth rates were observed at higher temperature. High variability of strain-specific responses has probably masked the overall effect of the treatments, highlighting the importance of testing the effect of environmental factors on many strains. It is likely that A. ostenfeldii populations can persist in the future, because suitable strains, which are able to germinate and grow well at potential future climate conditions, are part of the highly diverse cyst population. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://doi.org/10.5061/dryad.c8c83nr.

Consistency in microbiomes in cultures of Alexandrium species isolated from brackish and marine waters.

Phytoplankton and bacteria interactions have a significant role in aquatic ecosystem functioning. Associations can range from mutualistic to parasitic, shaping biogeochemical cycles and having a direct influence on phytoplankton growth. How variations in phenotype and sampling location, affect the phytoplankton microbiome is largely unknown. A high-resolution characterization of the bacterial community in cultures of the dinoflagellate Alexandrium was performed on strains isolated from different geographical locations and at varying anthropogenic impact levels. Microbiomes of Baltic Sea Alexandrium ostenfeldii isolates were dominated by Betaproteobacteria and were consistent over phenotypic and genotypic Alexandrium strain variation, resulting in identification of an A. ostenfeldii core microbiome. Comparisons with in situ bacterial communities showed that taxa found in this A. ostenfeldii core were specifically associated to dinoflagellate dynamics in the Baltic Sea. Microbiomes of Alexandrium tamarense and minutum, isolated from the Mediterranean Sea, differed from those of A. ostenfeldii in bacterial diversity and composition but displayed high consistency, and a core set of bacterial taxa was identified. This indicates that Alexandrium isolates with diverse phenotypes host predictable, species-specific, core microbiomes reflecting the abiotic conditions from which they were isolated. These findings enable in-depth studies of potential interactions occurring between Alexandrium and specific bacterial taxa.

Intraspecific trait variation and trade-offs within and across populations of a toxic dinoflagellate.

Intraspecific trait diversity can promote the success of a species, as complementarity of functional traits within populations may enhance its competitive success and facilitates resilience to changing environmental conditions. Here, we experimentally determined the variation and relationships between traits in 15 strains of the toxic dinoflagellate Alexandrium ostenfeldii derived from two populations. Measured traits included growth rate, cell size, elemental composition, nitrogen uptake kinetics, toxin production and allelochemical potency. Our results demonstrate substantial variation in all analysed traits both within and across populations, particularly in nitrogen affinity, which was even comparable to interspecific variation across phytoplankton species. We found distinct trade-offs between maximum nitrogen uptake rate and affinity, and between defensive and competitive traits. Furthermore, we identified differences in trait variation between the genetically similar populations. The observed high trait variation may facilitate development and resilience of harmful algal blooms under dynamic environmental conditions.

The potential for dispersal of microalgal resting cysts by migratory birds.

Most microalgal species are geographically widespread, but little is known about how they are dispersed. One potential mechanism for long-distance dispersal is through birds, which may transport cells internally (endozoochory) and deposit them during, or in-between, their migratory stopovers. We hypothesize that dinoflagellates, in particular resting stages, can tolerate bird digestion; that bird temperature, acidity, and retention time negatively affect dinoflagellate viability; and that recovered cysts can germinate after passage through the birds' gut, contributing to species-specific dispersal of the dinoflagellates across scales. Tolerance of two dinoflagellate species (Peridiniopsis borgei, a warm-water species and Apocalathium malmogiense, a cold-water species) to Mallard gut passage was investigated using in vitro experiments simulating the gizzard and caeca conditions. The effect of in vitro digestion and retention time on cell integrity, cell viability, and germination capacity of the dinoflagellate species was examined targeting both their vegetative and resting stages. Resting stages (cysts) of both species were able to survive simulated bird gut passage, even if their survival rate and germination were negatively affected by exposure to acidic condition and bird internal temperature. Cysts of A. malmogiense were more sensitive than P. borgei to treatments and to the presence of digestive enzymes. Vegetative cells did not survive conditions of bird internal temperature and formed pellicle cysts when exposed to gizzard-like acid conditions. We show that dinoflagellate resting cysts serve as dispersal propagules through migratory birds. Assuming a retention time of viable cysts of 2-12 h to duck stomach conditions, cysts could be dispersed 150-800 km and beyond.

A planktonic diatom displays genetic structure over small spatial scales.

Marine planktonic microalgae have potentially global dispersal, yet reduced gene flow has been confirmed repeatedly for several species. Over larger distances (>200 km) geographic isolation and restricted oceanographic connectivity have been recognized as instrumental in driving population divergence. Here we investigated whether similar patterns, that is, structured populations governed by geographic isolation and/or oceanographic connectivity, can be observed at smaller (6-152 km) geographic scales. To test this we established 425 clonal cultures of the planktonic diatom Skeletonema marinoi collected from 11 locations in the Archipelago Sea (northern Baltic Sea). The region is characterized by a complex topography, entailing several mixing regions of which four were included in the sampling area. Using eight microsatellite markers and conventional F-statistics, significant genetic differentiation was observed between several sites. Moreover, Bayesian cluster analysis revealed the co-occurrence of two genetic groups spread throughout the area. However, geographic isolation and oceanographic connectivity could not explain the genetic patterns observed. Our data reveal hierarchical genetic structuring whereby despite high dispersal potential, significantly diverged populations have developed over small spatial scales. Our results suggest that biological characteristics and historical events may be more important in generating barriers to gene flow than physical barriers at small spatial scales.

Delineating closely related dinoflagellate lineages using phylotranscriptomics.

Recently radiated dinoflagellates Apocalathium aciculiferum (collected in Lake Erken, Sweden), Apocalathium malmogiense (Baltic Sea) and Apocalathium aff. malmogiense (Highway Lake, Antarctica) represent a lineage with an unresolved phylogeny. We determined their phylogenetic relationships using phylotranscriptomics based on 792 amino acid sequences. Our results showed that A. aciculiferum diverged from the other two closely related lineages, consistent with their different morphologies in cell size, relative cell length and presence of spines. We hypothesized that A. aff. malmogiense and A. malmogiense, which inhabit different hemispheres, are evolutionarily more closely related because they diverged from a marine common ancestor, adapting to a wide salinity range, while A. aciculiferum colonized a freshwater habitat, by acquiring adaptations to this environment, in particular, salinity intolerance. We show that phylotranscriptomics can resolve the phylogeny of recently diverged protists. This has broad relevance, given that many phytoplankton species are morphologically very similar, and single genes sometimes lack the information to determine species' relationships.

Sample Preservation, DNA or RNA Extraction and Data Analysis for High-Throughput Phytoplankton Community Sequencing.

Phytoplankton is the basis for aquatic food webs and mirrors the water quality. Conventionally, phytoplankton analysis has been done using time consuming and partly subjective microscopic observations, but next generation sequencing (NGS) technologies provide promising potential for rapid automated examination of environmental samples. Because many phytoplankton species have tough cell walls, methods for cell lysis and DNA or RNA isolation need to be efficient to allow unbiased nucleic acid retrieval. Here, we analyzed how two phytoplankton preservation methods, three commercial DNA extraction kits and their improvements, three RNA extraction methods, and two data analysis procedures affected the results of the NGS analysis. A mock community was pooled from phytoplankton species with variation in nucleus size and cell wall hardness. Although the study showed potential for studying Lugol-preserved sample collections, it demonstrated critical challenges in the DNA-based phytoplankton analysis in overall. The 18S rRNA gene sequencing output was highly affected by the variation in the rRNA gene copy numbers per cell, while sample preservation and nucleic acid extraction methods formed another source of variation. At the top, sequence-specific variation in the data quality introduced unexpected bioinformatics bias when the sliding-window method was used for the quality trimming of the Ion Torrent data. While DNA-based analyses did not correlate with biomasses or cell numbers of the mock community, rRNA-based analyses were less affected by different RNA extraction procedures and had better match with the biomasses, dry weight and carbon contents, and are therefore recommended for quantitative phytoplankton analyses.